Field-induced magnetic phases and electric polarization in LiNiP04 
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Neutron difTraction is used to probe the (_ff, T) phase diagram of magneto-electric (ME) LiNiP04 
for magnetic fields along the c-axis. At zero field the Ni spins order in two antiferromagnetic phases. 
One has commensurate (C) structures and general ordering vectors kc = (0,0,0), the other one 
is incommensurate (IC) with k/c = (0,q,0). At low temperatures the C order collapses above 
fioH = 12 T and adopts an IC structure with modulation vector parallel to k/c- We show that C 
order is required for the ME effect and establish how electric polarization results from a field-induced 
reduction of the total magneto-elastic energy. 

PACS numbers: 75.25.+Z, 75.30.Gw, 75.80.-|-q 



Materials with both magnetic and electric order as 
found in magneto-electric (ME) multiferroics have re- 
ceived growing interest in recent years [J, 0; Q • 1^ 
pected that the coupling of magnetic and electric order 
in multiferroics will be of technological use, but also lead 
to rich physics with multi-order phase transitions 0] and 
excitations such as electromagnons @. Often ferro- 
electric and magnetic phases have very different ordering 
temperatures, suggesting that they arc driven by differ- 
ent microscopic interactions, but for some they coincide 
and ferroelectricity is generated by magnetic long-range 
order In the lithium ortho-phosphates, LiMP04, 

M=Mn, Fe, Co, Ni, a strong ME effect is observed in the 
antiferromagnetic phases. Mercier [l^l explained the tem- 
perature dependence of the ME coefficients for LiMnP04, 
LiCoP04 and to lesser extent for LiFcP04 with a micro- 
scopic model adapted from Cr203 [llj. However, the ME 
effect in LiNiP04 differs from the other lithium phos- 
phates and was not modelled as succesfuUy. In the present 
study we determine the field induced magnetic structures 
in LiNiP04 for fields H| |c and correct the existing picture 
of the zero field structures. Combining symmetry argu- 
ments and microscopic calculations similar to H , l3l we 
use the detailed information of the magnetic structures 
to quantify the ME properties of LiNiP04 and show how 
magnetic fields may lead to electric polarization. 

In LiNiP04, an electric polarization along the a-axis 
(c-axis) is generated when a magnetic field is applied 
along the c-axis (a-axis). This occurs below T = 20.8 K, 
where at zero field the material undergoes a first-order 
transition from a low temperature commensurate (C) 
antiferromagnetic phase with a general ordering vector 
kc = (0,0,0) for each the four Ni-spin sublattices (cf. 
Fig. [1^) to an incommensurate (IC) phase with kjc = 
(0,g,0) and 0.07 < g < 0.155 For fields H||c, the 
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FIG. 1: (Color online) (a) Ni^^ spin configuration of the zero- 
field C structure. Ni-positions are labelled i according to their 
positions ri (i = 1...4) but shifted (—0.25, —0.25, 0) compared 
to the values given in the text, (b-c) Projected C structure 
at zero and finite field H||c seen along the &-axis. Spin angles 
are exaggerated for clarity. The applied field cause the Ni2+ 
moments to rotate resulting in a magnetization M of the crys- 
tal, as described in the text, (d) Linearly polarized (LP) IC 
magnetic structure at zero field seen along the a-axis. 



magnetization measurements provide evidence for several 
phase transitions between hqH ~ 12 and 22 T lj|. Also 
IC magnetic structures have been discussed as possible 
explanation of hysteresis observed in the ME coefficients 
at high magnetic fields along the a-axis [l^ [3l ■ 

To understand the ME effect in LiNiP04, we have 
studied the (_ff , T) phase diagram and the magnetic 
structures for fields H||c up to ^qH = 14.7 T. We 
first present the (i7, T) phase diagram and show that 
off-diagonal single-ion anisotropics and Dzyaloshinsky- 
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Moriya (DM) interactions allowed by symmetry are con- 
sistent with the observed magnetic structures and lead to 
staggered magnetic moments in applied magnetic fields. 
Then we establish that electric polarization is only al- 
lowed in the field-induced C structure, but not in the 
high temperature IC and the zero-field C structures. Fi- 
nally we show that the electric polarization is driven by 
the magnetic symmetry and propose a model accounting 
for the temperature dependence of the ME constants. 

Measurements were performed on a high-quality 0.4 g 
single crystal. Zero field measurements were performed 
in a closed cycle cryostat on a four circle goniometer 
at the Tries single crystal diffractometer, using neutron 
wavelength A = 1.18 A for the C structure determination 
at T = 5 K, and A = 2.318 A for the IC structure at 
T = 21 K. For diffraction measurements on the triple 
axis spectrometer RITA-II the sample was mounted in a 
15 T magnet with the vertical field along the crystallo- 
graphic c-axis. The {H, T) phase diagram was determined 
using neutrons with A = 4.04 A. The high field magnetic 
structure was studied with A = 2.02 A neutrons. 

Symmetry properties: LiNiP04 crystallizes in the or- 
thorhombic Pnma (No. 62) crystal structure with lat- 
ti ce p arameters a = 10.02 A, b ~ 5.83 A and c = 4.66 
A[l7|. The magnetic Ni^+ ions with spin 5 = 1 are situ- 
ated on 4(c) sites forming buckled planes perpendicular 
to the a-axis. The positions of the four Ni^+ in each unit 
cell are: ri = (0.275,0.25,0.98), ra = (0.775,0.25,0.52), 
ra = (0.725,0.75,0.02) and r4 = (0.225,0.75,0.48), as 
shown in Fig. [TJi. The low crystal-field symmetry in 
LiNiP04 leads to a magnetic susceptibility tensor that 
contains staggered off-diagonal terms, Xac and Xca- 

This 

allows for a single-ion anisotropy of the following type: 

nrani ( _L C- ^^\ 

and two DM interactions: 

n/-DM 7-» f QC na qc na qc na i oc Qa\ ^ ' 

Phase diagram: The (H, T) phase diagram for fields 
H| |c up to ^J.oH = 14.7 T is shown in Fig. H The C phase 
is characterized by commensurate Bragg peaks associated 
with ordering vector kc, such as (0,1,0) whose T and 
H dependence are shown in Fig. O^-c. A sudden drop 
in intensity of the (0, 1,0) peak indicates the collapse of 
the C phase in a first-order phase transition. Between 
T = 10 K and 18 K, the C phase extends to higher 
fields, leading to a dome-shaped (iJ, T) phase diagram. 
The C phase is enclosed by an IC phase with a magnetic 
ordering wave-vector kjc, appearing e.g. at (0, 1 -I- g,0). 
Fig. [3^-d shows the temperature dependence of q and the 
intensity of the IC (0, 1 -I- q, 0) peaks for different fields. 

Staggered crystal-fields and DM interactions: The zero 
field C structure belongs to a single irreducible repre- 
sentation of kc, determined from 112 magnetic peaks 



at 5 K. The magnetic moments are nearly parallel 
to the c-axis with mc = (0.3(1), 0, 2.2(2))/LtB. The c- 
component, m^, has a (+,+,—,—) order and the a- 
componcnt, to^, a (-f , — , — , -I-) order on the sites with 
increasing i = 1...4 (Fig. [T^-b). Earlier structural analy- 
sis using powder diffraction [11, [l^ found , but not 
the smaller m^. The presence of may be explained 
by single ion anisotropics and DM interactions: Insert- 
ing itLq of (+,+,—,—) symmetry into Eqs. [T] and [21 
we find that JC^S' = -D,.^S{S1 - - + 5|) and 
JfDM ^ 2^(55^ - + 5|), which both favor 

that TO^ is of (+,—,—,-!-) symmetry. 

Field-induced staggered moments: Magnetic fields 
along the c-axis induce antiferromagnetic (1, 1,0) peaks 
that grow as {fioH)"^ in the C phase (Fig. [3b). High- 
energy (100 keV) X-ray diffraction at the BW5 beam 
line at HASYLAB, DESY detected no field dependent 
signal at (1, 1,0), showing that the neutron signal is of 
magnetic origin. Structural refinements reveal that the 
(1,1,0) peak reffects an additional staggered magnetic 
a-component, with (-I-, — , -|-, — ) symmetry, and an 
ordered moment increasing linearly with field to a value 
of nist ~ 0.17/iB at 12 T. Assuming that, to first order, 
the field rotates the magnetic moments without changing 
their magnitude, the magnetic structure (Fig. [TJ;) has a 
total ferromagnetic moment of (0,0, 0.6m"() per unit 
cell and the magnetization is M « (0, 0, 3.2)G at 12 T. 
This is consistent with bulk measurements [20[ showing 
that the magnetization grows almost linearly with H and 
is ~0.03 G at 0.1 T. At 40 K, where the system is param- 
agnetic, there is still field induced (1,1,0) intensity (Fig. 
[5^). We interpret this as staggered magnetic fields at the 
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FIG. 2: (Color online) (a) {H, T) phase diagram of LiNiP04 
for H||c including a commensurate (C) and an incommensu- 
rate (IC) antiferromagnetic phase, and a paramagnetic (P) 
phase, which at zero field supports short range fluctuations 
up to r = 40 K [l3|. (b-c) T and H dependence of the 
(0,1,0) magnetic Bragg peak intensity at three temperatures 
and fields, respectively. 
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FIG. 3: (Color online) (a-b) Position and integrated intensity 
of the (0,1+5,0) IC peak as function of temperature for dif- 
ferent fields H||c. (c-d) Scattering intensity for wave- vectors 
(0, fc, 0) at ^0^^ = 13.5 T and 14.7 T. At 13.5 T, the IC peak 
disappears near T = 15 K where the system enters the C 
phase (cf. Fig. [2^). (e-f) Background subtracted peak intensi- 
ties of the (1,1,0) peak as function of field and temperature. 



Ni position due to ofF-diagonal elements of the local sus- 
ceptibility tensor or the DM interactions, as prev iously 
observed in antiferromagnetic S=^ chains 
The IC magnetic order at zero field 



22|. 
is a trans- 
versely polarized coUinear spin-density wave belong- 
ing to a single representation of kjc. The structure is 
shown in Fig. [TJi and consists of magnetic moments 
mic = (0.0(2), 0.0(1), 1.2(2))/.iB that are ordered with 
(-I-, -t-, — /?, — /3) symmetry, where /3 = e"*'^'' describes the 
IC modulation along b. The IC component of the high- 
field structure is found by structural analysis of 83 mag- 
netic peaks at 14.7 T and 2.3 K. The best fit results in a 
similar structure as at zero field, but with increased am- 
plitude mic = (0, 0, 2.0(2))/xb. However, the data does 
not exclude an elliptically polarized (EP) IC structure 
with (-I-, —j3, —(i) components along a and c. An EP 
IC structure at 14.7 T and 2.3 K is directly supported 
by analysis of higher order harmonics and indirectly by 
a mean field calculation predicting a phase boundary be- 
tween the EP IC structure and a high temperature lin- 
early polarized (LP) IC structure around 15 K at 14.7 
T [26[. Coexisting with the IC order is a field- induced 



(1,1,0) intensity (Fig. [3^-f) signaHng a (-1-,—,+,—) C- 
type moment along the a-axis of approximately 0.11/.ts. 

Phenomenology: The zero field C structure breaks in- 
version symmetry, but is invariant under 2^ (180° screw 
axis along b), thus preventing electric polarization per- 
pendicular to b. However, for H||c the invariance under 
2b is broken and electric polarization is allowed. The LP 
IC structures leave at least one point of inversion invari- 
ant and do not allow for electric polarization - even in the 
presence of the C staggered moments. This is consistent 
with a more formal treatment developed by Harris [23 1. 

Magneto- electric effect: The main features of the ME 
effect in LiNiP04 can be explained by connections be- 
tween super-exchange (SE), DM spin interactions and 
elastic distortions. At zero field |Si| = IS2I = IS3I ^ 
S4I = (5*), the thermal mean value of the spin opera- 
tor, and the angles between Si and S2, and between S3 
and S4, are identical 612 = ^34 = 9. In the C phase a 
magnetic field H||c rotates the spins as shown in Fig. [41;. 
Here 6*12 = 6* + M and 6*34 = 6* - A^, and A6' is pro- 
portional to the magnetization XcHz, if we assume that 
jSil = IS2I = IS3I = IS4I even at non-zero fields. The SE 
energy for Oiff^^ — J12S1 ■ S2 -I- J34S3 • S4 in this spin 
configuration is 



£?f.34 + J3,){S)'{1 - -{9' + A9^)) 



(J12- .hi){sfeM. 



(3) 



The SE energy ([Sj can be lowered by a uniform dis- 
placement of exchange mediating ions such as the trans- 
lation of all P04-tetrahcdra along a by a small distance 
X (Fig. [3)3) . The symmetry of the Ni-O-P-O-Ni exchange 
paths implies that a uniform translation of the tetrahe- 
dra, leading to an electric polarization Px along a, simul- 
taneously increases J12 and reduces J34, or vice versa. To 
first order J12 = J -\-5 and J34 = J — 5, where i5 = \xX for 
small values of x. Introducing an elastic energy exx'^ for 
the tetrahedra displacements gives a SE-elastic interac- 
tion energy —2\x{S)^9lS.9x-\-exx'^, which is minimum for 
X = Xx{Sf9A9ex^. Noting that P^ (x x and AO cx Xc^^z 
we obtain an electrical polarization Px oc {S)'^XcHz, 
and thereby a ME coefficient axz oc e~^{S)'^Xc- An equiv- 
alent expression for a^z can also be obtained from the 
DM interaction term Hf^. These ME coefficients are 
similar to the phenomenological expressions suggested by 
Rado [3] but are here established from a microscopical 
point of view related to [ll|. Figure [H; compares the tem- 
perature dependence of the measured ME coefficient axz 
[l3| to e~^{Syxc (dashed line), assuming a constant elas- 
tic coefficient ex and using the magnetic order parameter 
(S) determined in [l^ and the magnetic susceptibility 
Xc from [iO, [2^. A more elaborate calculation of axz, 
assuming the angle difference, A^, fixed at the low tem- 
perature value, while the spins have non-identical lengths 
at finite temperatures, gives a significantly better agree- 
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FIG. 4: (Color online) (a) Positions of Ni (large blue (dark) 
circles), O (small yellow (bright) circles) and the PO4- 
tetrahedra (triangles) in zero field, (b) Same as (a), but for 
H||c. P04-tetrahedra are assumed to shift downwards with 
X (arrows), giving an electric polarization Px and changing 
the SE interaction as explained in the text, (c) Measured (cir- 
cles) [13 and calculated ME coefficient a^z for Hllc, assuming 
identical spin lengths (dashed line) and different spin lengths 
(solid line), (d) Same as (c) for a,x with H||a. Insets in (c) 
and (d) show the assumed angles between Si, S2, S3 and S4. 

ment with the experimental data (soUd line in Fig. 2);). 
Here the expressions for SE and DM are not equivalent 
and both terms are needed in the best fit to the data. 

To explain the elastic distortions in the C phase for 
H||a we first assume, as for H||c, that the magnetization 
of the sample results from a rotation of the magnetic mo- 
ments. This way we obtain the C spin structure sketched 
in Fig. HJi. Using similar arguments, now on the pairs 
(81,84) and (82,83). we find for identical spin lengths 
a ME coefficient Uzx oc (I'^{S)^Xa (dashed fine), which 
is compared to the measured ME coefficient Uzx in 
Fig. mi. Once again the elaborate calculation (solid line) 
improves the agreement with the experimental data. 

The proposed mechanism for ME distortions is not ef- 
fective for H||b in the C phase, nor for any field direc- 
tion in the LP IC phase. In the former case, all spins 
will cant with the same amount in the field direction and 
have the same lengths, leading to no energy differences 
between any pair of spins and no magneto-electricity. In 
the LP IC case, the observed spin structures are superpo- 
sitions of C and IC components, which are uncoupled in 
the energy terms because of translational symmetry. The 
symmetries of the C and the IC components considered 
separately do not produce the needed energy differences 
and ME distortions are therefore not induced. 



Conclusions: The symmetries of the established mag- 
netic structures do not support an electric polarization in 
the C- and zero field IC phase. Applying a magnetic field 
along c in the C phase creates a polar axis and allows for 
electric polarization. Symmetry analysis show that elec- 
tric polarization is possible in the C phase structure, but 
not in the LP IC phase. A microscopic model explains the 
temperature dependence of the ME constants, providing 
evidence that the electric polarization in LiNiP04 results 
from field induced changes in the magnetic structure. 
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